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Abstract  Every  three  years  the  IAU/IAG  Working  Group  on  Cartographic  Coor¬ 
dinates  and  Rotational  Elements  revises  tables  giving  the  directions  of  the  poles  of 
rotation  and  the  prime  meridians  of  the  planets,  satellites,  minor  planets,  and  com¬ 
ets.  This  report  introduces  improved  values  for  the  pole  and  rotation  rate  of  Pluto, 
Charon,  and  Phoebe,  the  pole  of  Jupiter,  the  sizes  and  shapes  of  Saturn  satellites  and 
Charon,  and  the  poles,  rotation  rates,  and  sizes  of  some  minor  planets  and  comets.  A 
high  precision  realization  for  the  pole  and  rotation  rate  of  the  Moon  is  provided.  The 
expression  for  the  Sun’s  rotation  has  been  changed  to  be  consistent  with  the  planets 
and  to  account  for  light  travel  time 
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1  Introduction 

The  IAU  Working  Group  on  Cartographic  Coordinates  and  Rotational  Elements  of 
the  Planets  and  Satellites  was  established  as  a  consequence  of  resolutions  adopted 
by  Commissions  4  and  16  at  the  IAU  General  Assembly  at  Grenoble  in  1976.  The 
Working  Group  became  a  joint  working  group  of  the  IAU  and  the  International 
Association  of  Geodesy  (IAG)  in  1985.  Now  within  the  IAU  the  working  group  is 
a  joint  working  group  of  Divisions  I  and  III,  and  not  part  of  commissions.  The  first 
report  of  the  Working  Group  was  presented  to  the  General  Assembly  at  Montreal  in 
1979  and  published  in  the  Trans.  IAU  17B,  72-79,  1980.  The  report  with  appendices 
was  published  in  Celestial  Mechanics  22,  205-230,  1980.  The  guiding  principles  and 
conventions  that  were  adopted  by  the  Group  and  the  rationale  for  their  acceptance 
were  presented  in  that  report  and  its  appendices.  The  second  report  of  the  Work¬ 
ing  Group  was  published  in  the  Trans.  IAU  18B,  151-162,  1983,  and  also  in  Celestial 
Mechanics  29,  309-321, 1983.  In  2003  the  name  of  the  Working  Group  was  shortened 
to  the  Working  Group  on  Cartographic  Coordinates  and  Rotational  Elements.  The 
table  summarizes  the  references  to  all  the  reports. 


Report 

General  Assembly 

Celestial  Mechanics  and  Dynamical  Astronomy 

1 

Montreal  in  1979 

22, 205-230  (Davies  et  al.  1980). 

2 

Patras  in  1982 

29,  309-321  (Davies  et  al.  1983). 

3 

New  Delhi  in  1985 

39, 103-113  (Davies  et  al.  1986). 

4 

Baltimore  in  1988 

46, 187-204  (Davies  et  al.  1989). 

5 

Buenos  Aires  in  1991 

53,  377-397  (Davies  et  al.  1992). 

6 

Hague  in  1994 

63, 127-148  (Davies  et  al.  1996). 

7 

Kyoto  in  1997 

No  report 

8 

Manchester  in  2000 

82,  83-110  (Seidelmann  et  al.  2002). 

9 

Sydney  in  2003 

91,  203-215  (Seidelmann  et  al.  2005). 

10 

Prague  in  2006 

This  paper 

G.  A.  Krasinsky 

Institute  for  Applied  Astronomy,  St.  Petersburg,  Russia 
G.  Neumann 

NASA  Goddard  Space  Flight  Center,  Greenbelt,  MD,  USA 
J.  Oberst 

DLR  Berlin  Adlershof,  Berlin,  Germany 
P.  Stooke 

University  of  Western  Ontario,  London,  Canada 
E.  F.  Tedesco 

University  of  New  Hampshire,  Durham,  NH,  USA 
D.  J.  Tholen 

University  of  Hawaii,  Honolulu,  HI,  USA 
P  C.  Thomas 

Cornell  University,  Ithaca,  NY,  USA 
I.  P  Williams 

Queen  Mary,  University  of  London,  London,  UK 
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Reprints  and  preprints  of  the  previous  and  this  report  can  be  found  at  the  working 
group  web  site:  http://astrogeology.usgs.gov/Projects/WGCCRE.  Previous  reports  are 
also  available  at  the  web  site:  http://www.springerlink.com/content/100246. 

The  previous  report  introduced  and  recommended  a  consistent  system  of  coordi¬ 
nates  for  both  minor  planets  and  comets.  This  system  is  not  the  same  as  the  system 
for  planets  and  satellites,  which  is  not  being  changed.  Pluto  is  included,  as  in  the  past, 
in  the  system  of  planets.  It  is  recognized  that  the  existence  of  two  different  systems 
has  the  potential  for  confusion,  but  the  methods  required  for  minor  planets  and  com¬ 
ets  differ  sufficiently  to  justify  the  use  of  two  different  systems.  This  report  includes 
descriptions  of  the  two  systems;  one  for  planets  and  satellites  and  another  for  minor 
planets  and  comets.  The  use  of  a  uniform  system  for  minor  planets  and  comets  is 
highly  recommended. 

The  IAU  passed  Resolution  5A  at  the  General  Assembly  in  Prague  on  August  24, 
2006,  adopting  a  definition  of  a  planet,  which  changes  the  classification  of  Pluto  and 
some  other  solar  system  bodies.  This  report,  which  is  based  on  the  progress  of  the 
past  triennium,  has  retained  the  previous  classifications  of  solar  system  bodies.  Future 
versions  of  the  report  will  incorporate  changes  as  appropriate. 

At  the  IAU  General  Assembly  in  Prague  Brent  Archinal  was  elected  as  the  new 
and  third  chairman  of  this  Working  Group. 


2  Definition  of  rotational  elements  for  planets  and  satellites 

Planetary  coordinate  systems  are  defined  relative  to  their  mean  axis  of  rotation  and 
various  definitions  of  longitude  depending  on  the  body.  The  longitude  systems  of 
most  of  those  bodies  with  observable  rigid  surfaces  have  been  defined  by  refer¬ 
ences  to  a  surface  feature  such  as  a  crater.  Approximate  expressions  for  these  rota¬ 
tional  elements  with  respect  to  the  International  Celestial  Reference  Frame  (ICRF) 
(Ma  et  al.  1998)  have  been  derived.  The  ICRF  is  the  reference  frame  of  the  Inter¬ 
national  Celestial  Reference  System  and  is  itself  epochless.  There  is  a  small  (well 
under  0.1  arcsecond)  rotation  between  the  ICRF  and  the  mean  dynamical  frame 
of  J2000.0.  The  epoch  J2000.0,  which  is  JD  2451545.0  (2000  January  1  1200  hours), 
Barycentric  Dynamical  Time  (TDB),  is  the  epoch  for  variable  quantities,  which  are 
expressed  in  units  of  days  (86400  SI  seconds)  or  Julian  centuries  of  36525  days.  TDB 
is  the  reference  time  scale  for  time  dependent  variables.  TDB  was  clarified  in  defi¬ 
nition  at  the  IAU  General  Assembly  of  2006  in  Prague.  TDB,  sometimes  called 
Teph,  is  roughly  equivalent  to  Terrestrial  Time  (TT)  in  epoch  and  rate.  UTC,  TCB, 
and  TCG  differ  from  TT  in  epoch  and  rate.  For  more  information  on  reference 
systems  and  time  scales  see  Kovalevsky  and  Seidelmann  (2004),  http://www.iers.org, 
http://rorf.usno.navy.mil/ICRF/,  or  http://aa.usno.navy.mil/faq/docs/ICRS_doc.php. 

The  north  pole  is  that  pole  of  rotation  that  lies  on  the  north  side  of  the  invariable 
plane  of  the  solar  system.  The  direction  of  the  north  pole  is  specified  by  the  value  of  its 
right  ascension  ao  and  declination  So.  With  the  pole  so  specified,  the  two  intersection 
points  of  the  body’s  equator  and  the  ICRF  equator  are  ao  ±90°.  We  chose  one  of  these, 
ao  +  90°,  and  define  it  as  the  node  Q.  Suppose  the  prime  meridian  has  been  chosen  so 
that  it  crosses  the  body’s  equator  at  point  B.  We  then  specify  the  location  of  the  prime 
meridian  by  providing  a  value  for  W,  the  angle  measured  easterly  along  the  body’s 
equator  between  the  node  Q  and  the  point  B  (see  Fig.  1).  The  right  ascension  of  the 
point  Q  is  90°  +  ao  and  the  inclination  of  the  planet’s  equator  to  the  celestial  equator 
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Fig.  1  Reference  system  used 
to  define  orientation  of  the 
planets  and  their  satellites 


is  90°  —  ($o.  As  long  as  the  planet,  and  hence  its  prime  meridian,  rotates  uniformly, 
W  varies  nearly  linearly  with  time.  In  addition,  ao,  <$o,  and  W  may  vary  with  time  due 
to  a  precession  of  the  axis  of  rotation  of  the  planet  (or  satellite).  If  W  increases  with 
time,  the  planet  has  a  direct  (or  prograde)  rotation,  and,  if  W  decreases  with  time,  the 
rotation  is  said  to  be  retrograde. 

In  the  absence  of  other  information,  the  axis  of  rotation  is  assumed  to  be  normal 
to  the  mean  orbital  plane  of  the  planet  or  the  satellite;  Mercury1  and  most  of  the 
satellites  are  in  this  category.  For  many  of  the  satellites,  it  is  assumed  that  the  rotation 
rate  is  equal  to  the  mean  orbital  period  (i.e.  synchronous  rotation,  but  in  some  cases 
such  an  assumption  still  needs  to  be  validated). 

The  angle  W  specifies  the  ephemeris  position  of  the  prime  meridian.  For  planets 
or  satellites  without  any  accurately  observable  fixed  surface  features,  the  adopted 
expression  for  W  defines  the  prime  meridian  and  is  not  subject  to  correction  for  this 
reason.  However,  the  rotation  rate  may  be  redefined  for  other  reasons.  Where  possi¬ 
ble,  however,  the  cartographic  position  of  the  prime  meridian  is  defined  by  a  suitable 
observable  feature,  and  so  the  constants  in  the  expression  W  =  Wo  +  W  d,  where 
d  is  the  interval  in  days  from  the  standard  epoch,  are  chosen  so  that  the  ephemeris 
position  follows  the  motion  of  the  cartographic  position  as  closely  as  possible;  in  these 
cases  the  expression  for  W  may  require  emendation  in  the  future. 

Recommended  values  of  the  constants  in  the  expressions  for  ao,  <$o,  and  W,  in  celes¬ 
tial  equatorial  coordinates,  are  given  for  the  planets  and  satellites  in  Tables  1  and  2.  In 
general,  these  expressions  should  be  accurate  to  one-tenth  of  a  degree;  however,  two 
decimal  places  are  given  to  assure  consistency  when  changing  coordinates  systems. 
Zeros  have  sometimes  been  added  to  rate  values  (W)  for  computational  consistency 
and  are  not  an  indication  of  significant  accuracy.  Additional  decimal  places  are  given 
in  the  expressions  for  Mercury,  the  Moon,  Mars,  Saturn,  and  Uranus,  reflecting  the 
greater  confidence  in  their  accuracy.  Expressions  for  the  Sun  and  Earth  are  given 
to  a  similar  precision  as  those  of  the  other  bodies  of  the  solar  system  and  are  for 
comparative  purposes  only. 


1  Based  on  dynamical  arguments,  e.g.  Peale  (2006),  Mercury’s  obliquity  is  thought  to  be  no  more  than 
a  few  arcseconds;  however,  the  current  uncertainty  in  its  obliquity  is  several  degrees  (J.  L.  Margot 
2003,  private  communication). 
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Table  1  Recommended  values  for  the  direction  of  the  north  pole  of  rotation  and  the  prime  meridian 
of  the  Sun  and  planets 


uq,  <$o  are  ICRF  equatorial  coordinates  at  epoch  J2000.0. 

Approximate  coordinates  of  the  north  pole  of  the  invariable  plane  are  ao  =  273°  .85,  <5q  =  66°. 99. 
T  =  interval  in  Julian  centuries  (of  36525  days)  from  the  standard  epoch 
d  =  interval  in  days  from  the  standard  epoch. 

The  standard  epoch  is  JD  2451545.0,  i.e.  2000  January  1  12  hours  TDB 


Sun 

cq)  =  286°. 13 

Sq  =  63°  .87 

W  =  84°. 176  +  14°.184400(M 

(a) 

Mercury 

c*o  =  281.01  -  0.033T 

S0  =  61.45  -  0.005  T 

W  =  329.548  +  6.1385025d 

(b) 

Venus 

c*o  =  272.76 

S0  =  67.16 

W  =  160.20  -  1.4813688d 

(c) 

Earth 

cxq  =  0.00  -  0.641  T 

So  =  90.00  -  0.557T 

W  =  190.147  +  360.9856235 d 

Mars 

cxq  =  317.68143  -  0.1061  T 

So  =  52.88650  -  0.0609T 

W  =  176.630  +  350.89198226^ 

(d) 

Jupiter 

cxo  =  268.056595  -  0.006499 T  +  0°  .000117  sin  Ja  +  0°  .000938  sin  Jb 
+  0.001432  sin  Jc  +  0.000030  sin  Jd  +  0.002150  sin  Je 

So  =  64.495303  +  0.002413  T  +  0.000050  cos  Ja  +  0.000404  cos  Jb 
+  0.000617  cos  Jc  —  0.000013  cos  Jd  +  0.000926  cos  Je 

W  =  284.95  +  870.5366420^ 

(e) 

where  Ja  =  99° .360714  +  4850°.4046T,  Jb  =  175°.895369  +  1191°.9605T, 

Jc  =  300°. 323162  +  262°.5475T,  Jd  =  114°.012305  +  6070°.2476T,  Je  =  49°.511251  +  64°.3000T 


Saturn 

c*o  =  40.589  -  0.036T 

So  =  83.537  -  0.004T 

W  =  38.90  +  810.7939024^ 

(e) 

Uranus 

c*o  =  257.311 

S0  =  -15.175 

W  =  203.81  -  501.1600928^ 

(e) 

Neptune 

(xq  =  299.36  +  0.70  sin  N 
<50  =  43.46  -  0.51  cos  N 

W  =  253.18  +  536.3128492^  -  0.48  sin  N 

N  =  357.85  +  52.316T 

(e) 

Pluto 

ao  =  312.993 

So  =  6.163 

W  =  237.305  -  56.3625225 d 

(f) 

(a)  The  equation  W  for  the  Sun  is  now  corrected  for  light  travel  time  and  removing  the  aberration 
correction.  See  Appendix 

(b)  The  20°  meridian  is  defined  by  the  crater  Hun  Kal 

(c)  The  0°  meridian  is  defined  by  the  central  peak  in  the  crater  Ariadne 

(d)  The  0°  meridian  is  defined  by  the  crater  Airy-0 

(e)  The  equations  for  W  for  Jupiter,  Saturn,  Uranus  and  Neptune  refer  to  the  rotation  of  their  magnetic 
fields  (System  III).  On  Jupiter,  System  I  (Wj  =  67°. 1  +  877°. 900d)  refers  to  the  mean  atmospheric 
equatorial  rotation;  System  II  (Wjj  =  43° .3  +  870°.270d)  refers  to  the  mean  atmospheric  rotation 
north  of  the  south  component  of  the  north  equatorial  belt,  and  south  of  the  north  component  of  the 
south  equatorial  belt 

(f)  The  0°  meridian  is  defined  as  the  mean  sub-Charon  meridian 
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3  Coordinate  system  for  the  moon 

The  recommended  coordinate  system  for  the  Moon  is  the  mean  Earth/polar  axis  (ME) 
system.  This  is  in  contrast  to  the  principal  axis  (PA)  system,  sometimes  called  the  axis 
of  figure  system.  The  ME  system,  sometimes  called  the  mean  Earth/rotation  axes  sys¬ 
tem,  is  recommended  because  nearly  all  cartographic  products  of  the  past  and  present 
have  been  aligned  to  it  (Davies  and  Colvin  2000).  The  difference  in  the  coordinates 
of  a  point  on  the  surface  of  the  Moon  between  these  systems  is  approximately  860  m. 
In  past  reports  the  rotation  and  pole  position  for  the  Moon  have  been  given  for  the 
ME  system  using  closed  formulae.  For  convenience  for  many  users,  those  formulae 
are  repeated  here  in  Table  2.  However,  users  should  note  that  these  are  valid  only 
to  the  approximately  150  m  level  of  accuracy,  as  shown  e.g.  by  Konopliv  etal.  (2001, 
Fig.  3).  For  high  precision  work  involving  e.g.  spacecraft  operations,  high-resolution 
mapping,  and  gravity  field  determination,  it  is  recommended  that  a  lunar  ephemeris 
be  used  to  obtain  the  libration  angles  for  the  Moon  from  which  the  pole  position  and 
rotation  can  be  derived. 

Specifically,  the  NASA/JPL  DE403/LE403  (Developmental  Ephemeris  403/Lunar 
Ephemeris  403),  commonly  known  as  DE403,  is  considered  the  best  currently  avail¬ 
able  lunar  ephemeris.  See  the  website  http://ssd.jpl.nasa.gov/iau-comm4/  for  more 
information  on  the  DE403  and  how  to  obtain  a  copy.  The  development  of  a  new  JPL 
lunar  ephemeris  is  under  consideration  (E.  M.  Standish  et  al.  2007,  private  commun- 
ciation)  and,  if  it  does  become  available,  it  might  be  used  for  the  highest  possible 
accuracy.  Polynomial  representations  of  the  (Euler)  lunar  libration  angles  and  their 
rates  in  the  PA  system  are  stored  in  the  ephemeris  file.  These  three  libration  angles 
are: 

(a)  <p,  the  angle  along  the  ICRF  equator,  from  the  ICRF  X-axis  to  the  ascending 
node  of  the  lunar  equator; 

(b)  0 ,  the  inclination  of  the  lunar  equator  to  the  ICRF  equator;  and 

(c)  i/s,  the  angle  along  the  lunar  equator  from  the  node  to  the  lunar  prime  meridian. 

Coordinates  or  Euler  angles  in  the  ME  system  (vector  M)  can  be  rotated  to  the  PA 
system  (vector  P)  using  the  following  expression  (Konopliv  et  al.  2001,  p.  7): 

P  =  Rz  (63// . 8986)  Ry  (7  9" . 07 68)  Rx  (0" . 1462)  M  (1) 

Conversely,  coordinates  or  Euler  angles  in  the  PA  system  can  be  rotated  into  the  ME 
system  with: 

M  =  Rx(— 0//.1462)Ry(— 79//.0768)Rz(— 63//.8986)P  (2) 

where  Rx,Ry,  and  Rz  are  the  standard  rotation  matrices  for  right-handed  rotations 
around  the  X,  Y,  and  Z  axes  respectively. 

Therefore,  for  a  given  epoch,  the  user  should  obtain  (p,0 ,  and  ^  from  the  ephemeris 
file  and  store  them  as  the  vector  P,  apply  the  transformation  in  Eq.  2,  and  extract  the 
angles,  now  in  the  ME  system  from  the  vector  M.  These  angles  can  then  be  converted 
with: 


ao  =  <P  ~  90° 
<50  =  90°  -  e 
W= f 


4?)  Springer 
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giving  the  lunar  rotation  angles  in  the  standard  ao,  and  W  formulation  (of  Table  2) 

and  in  the  ME  system. 

Alternatively,  if  the  user  has  coordinates  for  a  point  in  ICRF  coordinates  (vector  I) 
that  they  wish  to  convert  to  ME  coordinates,  for  a  given  epoch  the  user  should  obtain 
(p,0,  and  x/r  from  the  ephemeris  file,  and  then  do  the  conversion: 

M  =  Rx(-0//.1462)Ry(-79//.0768)Rz(-63//.8986)Rz(^)Rx(<9)Rz(^)I  (3) 

with  M  now  being  the  coordinates  of  the  point  in  the  ME  system.  The  user  should  note 
that  the  numerical  values  for  the  rotations  in  Eqs.  1,  2,  and  3  are  specific  to  DE403 
and  are  different  for  past  and  future  ephemerides. 

Note  that  the  NASA/JPL  Navigation  and  Ancillary  Information  Facility  (NAIF) 
provides  software  and  files  to  facilitate  the  above  transformations.  This  includes  a  Plan¬ 
etary  Constants  Kernel  (PCK)  made  using  the  lunar  libration  information  extracted 
from  the  DE/LE  403  ephemeris,  and  a  special  lunar  frames  kernel  (FK)  providing  the 
specifications  and  data  needed  to  construct  the  PA  to  ME  system  transformation.  A 
new  version  of  the  PCK  will  also  be  provided  when  a  new  JPL  ephemeris  is  released. 
See  http://naif.jpl.nasa.gov  for  further  information.  Roncoli  (2005)  also  provides  useful 
information  on  lunar  constants  and  coordinates,  including  on  the  differences  between 
the  ME  and  PA  systems  and  on  the  DE403  ephemeris. 


4  Rotation  elements  for  planets  and  satellites 

The  rotation  rate  of  Saturn,  which  is  given  in  Table  1  is  based  on  Voyager  observations 
of  kilometer  wavelength  radio  signals.  Recent  Cassini  observations  (Giampieri  et  al. 
2006)  of  a  signal  in  Saturn’s  magnetic  field  gives  a  period  of  about  10  h  and  47  min, 
about  8  min  longer  than  the  previously  determined  period.  At  this  time  it  is  uncertain 
whether  this  is  the  true  rotation  rate  or  what  is  the  physical  mechanism  causing  the 
different  signals  (Stevenson  2006).  Hence,  the  rotation  rate  of  Saturn  has  not  been 
changed,  while  more  results  from  the  Cassini  mission  are  anticipated. 

The  rotation  rates  of  Uranus  and  Neptune  were  determined  from  the  Voyager 
mission  in  1986  and  1989.  The  uncertainty  of  those  rotation  rates  are  such  that  the 
uncertainty  of  the  actual  rotation  position  is  more  than  a  complete  rotation  in  each 
case. 

A  new  model  for  the  pole  position  and  rotation  of  Mars  has  been  proposed  by 
Konopliv  et  al.  (2006)  based  on  the  most  recent  spacecraft  data.  At  this  time,  fol¬ 
lowing  the  advice  of  the  NASA  Mars  Geodesy  and  Cartography  Working  Group 
(Duxbury  2006),  the  use  of  this  new  model  is  not  recommended  for  cartographic 
purposes.  This  is  for  a  number  of  reasons  including  that  for  the  immediate  future  the 
new  model  would  have  little  if  any  effect  on  cartographic  products,  and  also  that  it 
is  expected  to  be  significantly  changed  in  the  next  few  years  as  new  data  becomes 
available.  However,  users  with  high  accuracy  requirements,  such  as  Mars  gravity  field 
determination,  may  wish  to  consider  using  it. 

The  topographic  reference  surface  of  Mars  is  that  specified  in  the  final  MOLA 
Mission  Experiment  Gridded  Data  Record  (MEGDR)  Products  (Smith  etal.  2003). 
In  particular,  the  128  pixels/0  resolution,  radius  and  topographic  surfaces  are  recom¬ 
mended,  although  the  lower  resolution  versions  may  be  used  where  appropriate  and 
documented,  and  for  the  areas  poleward  of  ±88°  latitude. 
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For  Mercury  the  use  of  a  planetocentric,  east-positive  (right-handed)  system  was 
adopted  by  the  MESSENGER  project  more  than  6  years  ago  to  facilitate  geodetic 
analysis,  particularly  topography  and  gravity,  as  well  as  all  cartography.  The  Mariner 
10  mission  used  the  IAU/IAG  standard  system.  There  are  standard  transformations 
between  the  two  coordinate  sets.  For  the  Mars  Global  Surveyor  mission,  an  areocen- 
tric,  east-positive  system  was  used  despite  years  of  Mariner  4,  6,  7,  and  9  and  Viking 
data  mapped  with  the  IAU/IAG  standard  system. 


5  Rotational  elements  for  minor  planets  and  comets 

For  planets  and  satellites,  the  IAU  definition  of  north  pole  is  the  pole  that  lies  above 
the  invariant  plane  of  the  solar  system,  and  the  rotation  can  be  either  direct  or  ret¬ 
rograde.  For  minor  planets  and  comets,  given  substantial  indirect  evidence  for  large 
precession  of  the  rotational  poles  of  some  comets,  this  first  definition  needs  to  be 
rethought.  In  particular,  situations  exist  in  which  the  pole  that  is  clearly  “north”  in  the 
IAU  sense  precesses  over  several  decades  to  become  clearly  “south”  in  the  IAU  sense. 
Comet  2P/Encke,  which  is  likely  to  be  visited  by  spacecraft  in  the  foreseeable  future, 
is  a  prime  example  of  a  comet  for  which  very  large  precession  has  been  inferred. 

There  is  also  clear  evidence  for  excited  state  rotation  at  least  for  comet  lP/Ftalley 
and  minor  planet  Toutatis.  In  this  case,  the  angular  momentum  vector  moves  around 
on  the  surface  of  the  body.  The  rotational  spin  vector  describes  substantial  excursions 
from  the  angular  momentum  vector  during  the  course  of  the  7-day  periodicity  that  is 
seen  in  the  light  curve.  We  can,  therefore,  anticipate  cases  in  which  the  rotational  spin 
pole  moves  back  and  forth  between  north  and  south  on  a  time  scale  of  days.  Thus, 
there  is  the  issue  of  needing  to  change  our  definition  of  the  rotational  pole. 

The  choice  of  a  rotational  pole  for  a  body  in  simple  rotation  with  slow  precession 
is  straightforward.  One  can  choose  the  pole  that  follows  either  the  right-hand  rule  or 
the  left-hand  rule,  and  the  right-hand  rule  is  chosen  here.  This  would  be  the  “positive” 
pole  to  avoid  confusion  with  the  north-south  terminology.  Ideally  one  would  like  to 
choose  a  pole  for  excited  state  rotation  that  reduces  to  this  definition  as  the  rotational 
energy  relaxes  to  the  ground  state.  For  SAM  (short-axis  mode)  rotational  states,  it 
is  possible  to  define  a  body-fixed  axis  that  circulates  in  a  generally  complex  pattern 
about  the  angular  momentum  vector  and  this  approaches  the  simple  right-hand  rule 
definition  as  the  rotational  energy  relaxes  to  the  ground  state  of  simple  rotation. 
Presumably  the  appropriate  body-fixed  pole  is  the  axis  of  maximum  moment  of  iner¬ 
tia.  However,  the  definition  for  a  body  in  a  LAM  (long-axis  mode)  rotational  state  is 
not  so  obvious,  since  there  is  then  complete  rotation  about  the  long  axis  of  the  body 
as  well  as  rotation  about  a  short  axis.  In  this  case,  the  pole  should  be  taken  as  the 
minimum  moment  of  inertia  (the  long  axis  of  an  ellipsoid)  according  to  the  right-hand 
rule. 

As  specified  in  Sect.  6,  for  planets  and  satellites,  longitude  should  increase  mono- 
tonically  for  an  observer  fixed  in  inertial  space.  For  minor  planets  and  comets,  however, 
with  the  above  rule  for  poles,  this  definition  corresponds  always  to  a  left-hand  rule  for 
increasing  longitude,  since  the  concept  of  retrograde  rotation  is  no  longer  relevant. 
Therefore,  for  minor  planets  and  comets,  to  be  consistent  with  the  above  pole  defini¬ 
tion,  increasing  longitude  should  always  follow  the  right-hand  rule.  This  definition  is 
consistent  with  the  sense  of  increasing  longitude  used  for  Eros  by  Miller  et  al.  (2002), 
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but  is  inconsistent  with  the  sense  of  increasing  longitude  used  for  Eros  by  Thomas 
etal.  (2002). 

For  each  minor  planet  and  comet  the  positive  pole  of  rotation  is  selected  as  the 
maximum  or  minimum  moment  of  inertia  according  to  whether  there  is  short  or  long 
axis  rotational  state  and  according  to  the  right-hand  rule.  So  for  minor  planets  and 
comets  the  positive  pole  is  specified  by  the  value  of  its  right  ascension  ao  and  declina¬ 
tion  <$o.  With  the  pole  so  specified,  the  two  intersection  points  of  the  body’s  equator 
and  the  ICRF  equator  are  ao  ±  90°.  We  chose  one  of  these,  ao  +  90°,  and  define  it 
as  the  node  Q.  Suppose  the  prime  meridian  has  been  chosen  so  that  it  crosses  the 
body’s  equator  at  the  point  B.  We  then  specify  the  location  of  the  prime  meridian  by 
providing  a  value  for  W,  the  angle  measured  along  the  body’s  equator  between  the 
node  Q  and  the  point  B  in  a  right-hand  system  with  respect  to  the  body’s  positive 
pole  (see  Fig.  2).  The  right  ascension  of  the  point  Q  is  90°  +  ao  and  the  inclination 
of  the  body’s  equator  to  the  celestial  equator  is  90°  —  <$o.  As  long  as  the  planet,  and 
hence  its  prime  meridian,  rotates  uniformly,  W  varies  linearly  with  time  according  to 
the  right-hand  rule.  In  addition,  ao,  So,  and  W  may  vary  with  time  due  to  a  precession 
of  the  axis  of  rotation  of  the  body. 

The  angle  W  specifies  the  ephemeris  position  of  the  prime  meridian,  and  for 
minor  planets  or  comets  without  any  accurately  observable  fixed  surface  features, 
the  adopted  expression  for  W  defines  the  prime  meridian.  Where  possible,  however, 
the  cartographic  position  of  the  prime  meridian  is  defined  by  a  suitable  observable 
feature,  and  so  the  constants  in  the  expression  W  =  Wo  +  Wd,  where  d  is  the  interval 
in  days  from  the  standard  epoch,  are  chosen  so  that  the  ephemeris  position  follows  the 
motion  of  the  cartographic  position  as  closely  as  possible;  in  these  cases  the  expression 
for  W  may  require  emendation  in  the  future.  Table  3  gives  the  recommended  rota¬ 
tion  values  for  the  direction  of  the  positive  pole  of  rotation  and  the  prime  meridian 
of  selected  minor  planets  and  comets.  Values  are  given  for  objects  that  have  been 
imaged  by  spacecraft,  radar,  or  high  resolution  Earth  based  imaging  systems  with 
sufficient  resolution  to  establish  accurate  pole  orientation  and  rotation  rates.  Values 
are  not  given  for  objects  where  the  observations  are  limited  to  photometric  light 
curves. 


Fig.  2  Reference  system  used 
to  define  orientation  of  the 
minor  planets  and  comets 
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Table  3 

Recommended  rotation  values  for  the  direction  of  the  positive  pole  of  rotation  and  the 

prime  meridian  of  selected  minor  planets  and  comets 

d  is  the  interval  in  days  from  the  standard  epoch,  i.e.  J2000.0  = 
12  hours  TDB.  cxq,  <5q,  and  W  are  as  defined  in  the  text. 

JD  2451545.0,  i.e.  2000  January  1 

4 

Vesta 

cxq  =  301° 

«50  =  41° 

W  =  292°  +  1617°. 332776d 

(a) 

243 

Ida 

c*o  =  168°  .76 
$0  =  -2°  .88 

W  =  265°. 95  +  1864°. 6280070d 

(b) 

433 

Eros 

c*o  =  11°. 35  dfc  0.02 
<50  =  17°. 22  ±  0.02 

W  =  326°. 07  +  1639°. 38864745d 

(c) 

951 

Gaspra 

ao  =  9°. 47 
<$0  =  26°. 70 

W  =  83°  .67  +  1226°.9114850d 

(d) 

25143 

Itokawa 

aQ  =  90°. 53 
$0  =  -66°. 30 

W  =  000°  +  712°. 143d 

(e) 

9P/Tempel  1 

aQ  =  294° 

$0  =  73° 

W  =  252°. 63  +  212°.064<i 

(f) 

19P/Borrelly 

c*o  =  218°.5±3 

$0  =  —12°  .5  dr  3 
w  =  000°  +  390°. 0  d 

(e) 

(a)  The  0°  meridian  is  defined  by  the  Others  Regio  (informal  name) 

(b)  The  0°  meridian  is  defined  by  the  crater  Afon 

(c)  The  0°  meridian  is  defined  by  an  unnamed  crater 

(d)  The  0°  meridian  is  defined  by  the  crater  Charax 

(e)  Since  only  rotation  rate  information  is  available,  the  0°  meridian  is  currently  arbitrarily  defined 
with  Wo  =  0° 

(f)  The  0°  meridian  is  defined  by  a  350  m  diameter  unnamed  circular  feature  near  the  Deep  Impactor 
impact  site  (Thomas  et  al.  2007a) 


6  Definition  of  cartographic  coordinate  systems  for  planets  and  satellites 

In  mathematical  and  geodetic  terminology,  the  terms  ‘latitude’  and  ‘longitude’  refer 
to  a  right-hand  spherical  coordinate  system  in  which  latitude  is  defined  as  the  angle 
between  a  vector  passing  through  the  origin  of  the  spherical  coordinate  system  and 
the  equator,  and  longitude  is  the  angle  between  the  vector  and  the  plane  of  the 
prime  meridian  measured  in  an  eastern  direction.  This  coordinate  system,  together 
with  Cartesian  coordinates,  is  used  in  most  planetary  computations,  and  is  sometimes 
called  the  planetocentric  coordinate  system.  The  origin  is  the  center  of  mass. 

Because  of  astronomical  tradition,  planetographic  coordinates  (those  commonly 
used  on  maps)  may  or  may  not  be  identical  with  traditional  spherical  coordinates. 
Planetographic  coordinates  are  defined  by  guiding  principles  contained  in  a  resolu¬ 
tion  passed  at  the  fourteenth  General  Assembly  of  the  IAU  in  1970.  These  guiding 
principles  state  that: 

(1)  The  rotational  pole  of  a  planet  or  satellite  which  lies  on  the  north  side  of  the 
invariable  plane  will  be  called  north,  and  northern  latitudes  will  be  designated 
as  positive. 
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(2)  The  planetographic  longitude  of  the  central  meridian,  as  observed  from  a  direc¬ 
tion  fixed  with  respect  to  an  inertial  system,  will  increase  with  time.  The  range 
of  longitudes  shall  extend  from  0°  to  360°. 

Thus  west  longitudes  (i.e.,  longitudes  measured  positively  to  the  west)  will  be  used 
when  the  rotation  is  direct  and  east  longitudes  (i.e.,  longitudes  measured  positively 
to  the  east)  when  the  rotation  is  retrograde.  The  origin  is  the  center  of  mass.  Also 
because  of  tradition,  the  Earth,  Sun,  and  Moon  do  not  conform  to  this  definition. 
Their  rotations  are  direct  and  longitudes  run  both  east  and  west  180°,  or  east  360°. 

For  planets  and  satellites,  latitude  is  measured  north  and  south  of  the  equator; 
north  latitudes  are  designated  as  positive.  The  planetographic  latitude  of  a  point  on 
the  reference  surface  is  the  angle  between  the  equatorial  plane  and  the  normal  to  the 
reference  surface  at  the  point.  In  the  planetographic  system,  the  position  of  a  point 
(P)  not  on  the  reference  surface  is  specified  by  the  planetographic  latitude  of  the  point 
(. P' )  on  the  reference  surface  at  which  the  normal  passes  through  P  and  by  the  height 
(h)  of  P  above  P'. 

The  reference  surfaces  for  some  planets  (such  as  Earth  and  Mars)  are  ellipsoids  of 
revolution  for  which  the  radius  at  the  equator  (A)  is  larger  than  the  polar  semi-axis  (C). 

Calculations  of  the  hydrostatic  shapes  of  some  of  the  satellites  (Io,  Mimas,  Encela- 
dus,  and  Miranda)  indicate  that  their  reference  surfaces  should  be  triaxial  ellipsoids. 
Triaxial  ellipsoids  would  render  many  computations  more  complicated,  especially 
those  related  to  map  projections.  Many  projections  would  lose  their  elegant  and  pop¬ 
ular  properties.  For  this  reason  spherical  reference  surfaces  are  frequently  used  in 
mapping  programs. 

Many  small  bodies  of  the  solar  system  (satellites,  minor  planets,  and  comet  nuclei) 
have  very  irregular  shapes.  Sometimes  spherical  reference  surfaces  are  used  for  com¬ 
putational  convenience,  but  this  approach  does  not  preserve  the  area  or  shape  charac¬ 
teristics  of  common  map  projections.  Orthographic  projections  often  are  adopted  for 
cartographic  portrayal  as  these  preserve  the  irregular  appearance  of  the  body  without 
artificial  distortion.  A  more  detailed  discussion  of  cartographic  coordinate  systems 
for  small  bodies  is  given  in  Sect.  7  of  this  report. 

Table  4  gives  the  size  and  shape  parameters  for  the  planets.  In  that  table  average 
(AVG),  north  (N),  and  south  (S)  polar  radii  are  given  for  Mars.  For  the  purpose  of 
adopting  a  best-fitting  ellipsoid  for  Mars,  the  average  polar  radius  should  be  used— the 
other  values  are  for  comparison  only,  e.g.  to  illustrate  the  large  dichotomy  in  shape 
between  the  northern  and  southern  hemispheres  of  Mars.  In  applications  where  these 
differences  may  cause  problems,  the  earlier  recommended  topographic  shape  model 
for  Mars  should  probably  be  used  as  a  reference  surface. 

Table  5  gives  the  size  and  shape  of  satellites  where  known.  Only  brightnesses  are 
known  for  many  of  the  newly  discovered  satellites.  Poles  and  rotation  rates  are  also 
not  yet  known  for  the  new  discoveries,  so  those  satellites  are  not  listed. 

The  values  of  the  radii  and  axes  in  Tables  4  and  5  are  derived  by  various  methods 
and  do  not  always  refer  to  common  definitions.  Some  use  star  or  spacecraft  occulta- 
tion  measurements,  some  use  limb  fitting,  others  use  altimetry  measurements  from 
orbiting  spacecraft,  and  some  use  control  network  computations.  For  the  Earth,  the 
spheroid  refers  to  mean  sea  level,  clearly  a  very  different  definition  from  other  bodies 
in  the  Solar  System. 

The  uncertainties  in  the  values  for  the  radii  and  axes  in  Tables  4  and  5  are  generally 
those  of  the  authors,  and,  as  such,  frequently  have  different  meanings.  Sometimes 
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Table  4 

Size  and  shape  parameters  of  the  planets 

Planet 

Mean  Equatorial  Polar  radius 

radius  (km)  radius  (km)  (km) 

RMS 

deviation 

from 

spheroid 

(km) 

Maximum 

elevation 

(km) 

Maximum 

depression 

(km) 

Mercury 

2439.7  ±  1.0 

Same 

Same 

1 

4.6 

2.5 

Venus 

6051.8  ±  1.0 

Same 

Same 

1 

11 

2 

Earth 

6371.00  ± 

0.01 

6378.14  ±0.01 

6356.75  ±0.01 

3.57 

8.85 

11.52 

Mars 

3389.50  ± 

0.2 

3396.19  ± 

0.1 

AVG  3376.20 
±0.1 

N  3373.19  ± 

0.1 

S  3379.21  ± 

0.1 

3.0 

22.64  ± 

0.1 

7.55  ± 

0.1 

Jupiter* 

69911  ±  6 

71492  ± 4 

66854  ±  10 

62.1 

31 

102 

Saturn* 

58232  ±  6 

60268  ± 4 

54364  ±  10 

102.9 

8 

205 

Uranus* 

25362  ±  7 

25559  ± 4 

24973  ±  20 

16.8 

28 

0 

Neptune* 

Pluto 

24622  ±  19 

1195  +  5 

24764  ±  15 

Same 

24341  ±  30 

Same 

8 

14 

0 

*  The  radii  correspond  to  a  one-bar  surface 


they  are  standard  errors  of  a  particular  data  set,  sometimes  simply  an  estimate  or 
expression  of  confidence.  The  radii  and  axes  of  the  large  gaseous  planets,  Jupiter, 
Saturn,  Uranus,  and  Neptune  in  Table  4  refer  to  a  one-bar-pressure  surface.  The  radii 
given  in  the  tables  are  not  necessarily  the  appropriate  values  to  be  used  in  dynamical 
studies;  the  radius  actually  used  to  derive  a  value  of  J2  (for  example)  should  always 
be  used  in  conjunction  with  it. 

7  Cartographic  coordinates  for  minor  planets  and  comets 

For  large  bodies,  a  spherical  or  ellipsoidal  model  shape  has  traditionally  been  defined 
for  mapping,  as  in  our  past  reports.  For  irregularly  shaped  bodies  the  ellipsoid  is 
obviously  useless,  except  perhaps  for  dynamical  studies.  For  very  irregular  bodies, 
the  concept  of  a  reference  ellipsoid  ceases  to  be  useful  for  most  purposes.  For  these 
bodies,  topographic  shapes  are  usually  represented  by  a  grid  of  radii  to  the  surface  as 
a  function  of  planetocentric  latitude  and  longitude  (when  possible,  or  also  by  a  set  of 
vertices  and  polygons). 

Another  problem  with  small  bodies  is  that  two  coordinates  (i.e.  spherical  angular 
measures)  may  not  uniquely  identify  a  point  on  the  surface  of  the  body.  In  other 
words  it  is  possible  to  have  a  line  from  the  center  of  the  object  intersect  the  surface 
more  than  once.  This  can  happen  on  large  and  even  mostly  ellipsoidal  objects  such 
as  the  Earth,  because  of  such  features  as  overhanging  cliffs  and  natural  bridges  and 
arches.  However,  on  large  bodies  these  features  are  relatively  very  small  and  often 
ignored  at  the  scale  of  most  topographic  maps.  For  small  bodies  they  may  be  fairly 
large  relative  to  the  size  of  the  body.  Example  cases  are  on  Eros  (at  a  small  patch 
west  of  Psyche),  and  certainly  on  Kleopatra  (Ostro  2000),  possibly  on  Toutatis  near 
its  ‘neck’,  and  perhaps  near  the  south  pole  of  Ida,  some  radii  may  intersect  the  surface 
more  than  once.  Even  on  small  bodies  this  problem  is  usually  restricted  to  small  areas, 
but  it  still  may  make  a  planetocentric  coordinate  system  difficult  to  use.  Cartographers 
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Umbriel  584.7  ±  2.8  Same  Same 

Titania  788.9  ±  1.8  Same  Same 

Oberon  761.4  ±  2.6  Same  Same 
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always  have  ad  hoc  tricks  for  a  specific  map,  such  as  interpolating  across  the  problem 
area  from  areas,  which  are  uniquely  defined,  or  by  showing  overlapping  contours.  A 
Cartesian  or  other  coordinate  geometry  may  be  preferable  for  arbitrarily  complex 
shapes,  such  as  a  toroidal  comet  nucleus,  where  an  active  region  has  eaten  its  way 
through  the  nucleus.  Such  coordinate  geometries  may  also  be  useful  for  irregular 
bodies  imaged  only  on  one  side,  such  as  for  19P/Borrelly  and  81P/Wild  2. 

With  the  introduction  of  large  mass  storage  to  computer  systems,  digital  cartog¬ 
raphy  has  become  increasingly  popular.  Cartographic  databases  are  important  when 
considering  irregularly  shaped  bodies  and  other  bodies,  where  the  surface  can  be 
described  by  a  file  containing  the  coordinates  for  each  pixel.  In  this  case  the  reference 
sphere  has  shrunk  to  a  unit  sphere.  Other  parameters  such  as  brightness,  gravity,  etc., 
if  known,  can  be  associated  with  each  pixel.  With  proper  programming,  pictorial  and 
projected  views  of  the  body  can  then  be  displayed. 

Taking  all  of  this  into  account,  our  recommendation  is  that  longitudes  on  minor 
planets  and  comets  should  be  measured  positively  from  0  to  360  degrees  using  a  right- 
hand  system  from  a  designated  prime  meridian.  The  origin  is  the  center  of  mass,  to 
the  extent  known. 

Latitude  is  measured  positive  and  negative  from  the  equator;  latitudes  toward  the 
positive  pole  are  designated  as  positive.  For  regular  shaped  bodies  the  cartographic 
latitude  of  a  point  on  the  reference  surface  is  the  angle  between  the  equatorial  plane 
and  the  normal  to  the  reference  surface  at  the  point.  In  the  cartographic  system,  the 
position  of  a  point  (P)  not  on  the  reference  surface  is  specified  by  the  cartographic 
latitude  of  the  point  ( P' )  on  the  reference  surface  at  which  the  normal  passes  through 
P  and  by  the  height  ( h )  of  P  above  P'. 

For  irregular  bodies  orthographic  digital  projections  often  are  adopted  for  car¬ 
tographic  portrayal  as  these  preserve  the  irregular  appearance  of  the  body  without 
artificial  distortion.  These  projections  should  also  follow  the  right-hand  rule. 

Table  6  contains  data  on  the  size  and  shape  of  selected  minor  planets  and  comets. 
The  first  column  gives  the  effective  radius  of  the  body  and  an  estimate  of  the  accuracy 


Table  6  Size  and  shape  parameters  of  selected  minor  planets  and  comets 


Asteroid/  comet  Effective  radius  (km)  Radii  measured  along  principal  axes  (km) 


1  Ceres 

487.3  ±  1.8 

454.7  ±  1.6 

4  Vesta 

289  ±5 

280  ±5 

229  ±5 

243  Ida 

15.65  ±  0.6 

26.8 

12.0 

7.6 

253  Mathilde 

26.5  ±  1.3 

33 

24 

23 

433  Eros 

8.45  ±  0.02 

17.0 

5.5 

5.5 

951  Gaspra 

6.1  ±  0.4 

9.1 

5.2 

4.4 

4179  Toutatis 

2.13 

1.015 

0.85 

25143  Itokawa 

0.535 

0.294 

0.209 

lP/Halley 

8.0  ±  0.5 

4.0  ±  0.25 

4.0  ±  0.25 

9P/Tempel  1 

3.0  ±0.1 

3.7 

2.5 

19P/Borrelly 

4.22  ±  0.05 

3.5  ±  0.2 

- 

- 

81P/Wild  2 

1.975 

2.7 

1.9 

1.5 

(a)  An  oblate  spheroid 

(b)  The  maximum  and  minimum  radii  are  not  properly  the  values  of  the  principal  semi-axes,  they  are 
half  the  maximum  and  minimum  values  of  the  diameter.  Due  to  the  large  deviations  from  a  simple 
ellipsoid,  they  may  not  correspond  with  measurements  along  the  principal  axes,  or  be  orthogonal  to 
each  other 
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of  this  measurement.  This  effective  radius  is  for  a  sphere  of  equivalent  volume.  The 
next  three  columns  give  estimates  of  the  radii  measured  along  the  three  principal 
axes. 

The  uncertainties  in  the  values  for  the  radii  in  Table  6  are  generally  those  of  the 
authors,  and,  as  such,  frequently  have  different  meanings.  Sometimes  they  are  stan¬ 
dard  errors  of  a  particular  data  set,  sometimes  simply  an  estimate  or  expression  of 
confidence. 

The  radii  given  in  the  tables  are  not  necessarily  the  appropriate  values  to  be  used 
in  dynamical  studies;  the  radius  actually  used  to  derive  a  value  for  the  dynamical  form 
factor  (J2)  (for  example)  should  always  be  used  in  conjunction  with  it. 
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Appendix:  Changes  since  the  last  report 

This  appendix  summarizes  the  changes  that  have  been  made  to  the  tables  since  the 

2003  report  (Celestial  Mechanics  and  Dynamical  Astronomy  91, 203-215,  2005). 

1.  The  expression  for  the  Sun’s  rotation  has  been  changed  to  account  for  the  light 
travel  time  and  removing  the  aberration  correction.  The  value  in  W  of  84.176° 
has  replaced  84.10°.  The  value  84.10°  is  correct  for  the  case  where  d  is  meant  to 
be  TT,  when  the  light  arrives  at  the  Earth,  not  the  moment  when  the  light  left  the 
Sun.  The  84.176°  is  correct  for  the  time  when  the  light  left  the  Sun  without  the 
aberration  correction,  which  is  consistent  with  the  values  given  for  the  planets, 
whose  light  travel  time  is  not  as  constant  as  for  the  Sun. 

2.  The  pole  and  rotation  rate  of  Pluto  in  Table  1  and  of  Charon  in  Table  2  have  been 
improved  based  on  Tholen  and  Buie  1997.  The  pole  of  Jupiter  has  been  improved 
based  on  Jacobson  2005  private  communication  and  Jacobson  2002.  The  Mars 
value  of  Wo  has  not  changed,  but  the  correct  reference  is  given  here  as  Duxbury 
et  al.  2001. 

3.  The  expressions  of  the  pole  and  rotation  of  the  Moon  are  noted  in  Table  2  as 
being  for  low  precision  use.  An  algorithm  is  described  in  the  text  for  using  the 
JPL  DE403  lunar  ephemeris,  rotated  to  the  mean  Earth/polar  axis  system,  in 
order  to  obtain  the  pole  and  rotation  with  high  precision. 

4.  The  new  pole  and  rotation  models  of  Konopliv  et  al.  2006  for  Mars  and  Giampieri 
et  al.  2006  for  Saturn  are  noted  in  the  text,  but  not  recommended  for  general  use 
at  this  time. 

5.  The  pole  and  rotation  rate  of  Phoebe  have  been  updated  in  Table  2  and  the  size 
and  shape  have  been  updated  in  Table  5  (P.  C.  Thomas  private  communication, 
also  see  Porco  et  al.  2005). 

6.  The  pole  and  rotation  rate  of  Itokawa  has  been  added  to  Table  3  based  on  Fujiw- 
ara  et  al.  2006.  Values  for  9P/  Tempel  1  have  been  added  to  Tables  3  and  6  based 
on  Thomas  et  al.  2007a. 

7.  In  Table  5  the  sizes  and  shapes  of  Saturn  satellites  I,  II,  III,  IV,  and  V  have  been 
corrected  based  on  Thomas  et  al.  2007b.  The  mean  radius  of  I  Mimas  and  the  size 
and  shape  of  VIII  Iapeteus  are  from  P.  C.  Thomas  private  communication.  The 
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radius  of  Charon  has  been  updated  based  on  the  combination  of  Gulbis  et  al.  2006 
and  Sicardy  et  al.  2006. 

8.  In  Table  6  the  sizes  of  Ceres  and  Vesta  have  been  added  based  on  Thomas  et  al. 
2005  and  Thomas  et  al.  1997.  The  size  of  Itokawa  is  based  on  Fujiwara  et  al.  2006. 
Minor  planets  Mathilde,  Eros,  and  Toutatis,  and  comets  Halley,  Tempel  1,  and 
Wild  2  have  been  added.  Minor  planets  Kleopatra,  Golevka,  Nyx,  1998JM8,  and 
1998ML14  have  been  deleted,  because  they  were  modeled  from  low  resolution 
radar  data,  and  cannot  be  mapped  from  those  data. 

9.  The  time  epochs  have  been  restated  to  be  TDB,  rather  than  TT,  since  the  barycen- 
tric  time  is  technically  correct  for  these  tables.  TDB  and  TT  are  roughly  equivalent 
in  epoch  and  rate  and  to  the  accuracies  given  do  not  differ.  Also  the  JD  value  has 
been  given  before  the  calendar  date  for  clarity. 

10.  A  sign  error  has  been  corrected  in  the  equation  for  W  for  243  Ida  in  Table  3. 
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